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Starting from Feynman's Lagrangian description of quantum mechanics, we propose a method 
to construct explicitly the propagator for the Wigner distribution function of a single system. For 
general quadratic Lagrangians, only the classical phase space trajectory is found to contribute to 
the propagator. Inspired by Feynman's and Vernon's influence functional theory we extend the 
method to calculate the propagator for the reduced Wigner function of a system of interest coupled 
to an external system. Explicit expressions are obtained when the external system consists of a set 
of independent harmonic oscillators. As an example we calculate the propagator for the reduced 
Wigner function associated with the Caldeira-Legett model. 



I. INTRODUCTION 

To explain the dynamic behavior of a quantum mechanical system of interest one generally needs to extract the 
time dependence of the density matrix. Then, an arbitrary observable represented by a Hermitian operator can be 
obtained by taking the trace of the operator times the density matrix. 

In a classical description, however, the density matrix emerges as a continuous, time dependent phase space distri- 
bution function whereas physical observables are identified with ordinary functions. Multiplying the latter with the 
phase space distribution function and integrating over phase space yields the commonly known recipe to calculate all 
relevant expectation values. Both approaches can be linked if one addresses the Wigner function[l| which is nothing 
but the Wigner-Weyl transform of the quantum mechanical density matrix. Inheriting from the density matrix all 
necessary information to calculate any observables of interest, the Wigner satisfies the Wigncr-Liouville equation, the 
quantum mechanical analog of the classical Liouvillc equation. 

Alternative to the Hamiltonian description of quantum mechanics one may adopt Feynman's Lagrangian formu- 
lation of quantum mechanics @. This approach essentially relies on the the calculation of the Green's function or 
propagator appearing as a path integral, rather than on solving the dynamical equations equation for the density ma- 
trix. Moreover, also Feynman's approach exhibits a direct link to classical mechanics through the explicit occurrence 
of the classical action in the path integral. 

In this work, adopting exclusively Feynman's Lagrangian formulation, we construct the path integral for the prop- 
agator of the single particle Wigner function in section [TTJ In section Hill we extend the obtained result by computing 
also the path integral for the propagator of the reduced Wigner function of a system coupled to an external quantum 
system. This extension may be seen as the Wigner-Weyl formulation of the theory of influence functionals developed 
by Feynman and VernonQ- In section HVl we derive the explicit Wigner influence functional for a particle coupled 
to a set of independent harmonic oscillators. Finally, as a non trivial example, we calculate the Wigner function 
propagator of a particle described by the Caldeira-Legett model Q. 

II. SINGLE-PARTICLE PHASE SPACE PROPAGATOR 

The well-known Wigner distribution function [l| provides a phase space description of quantum mcchanics,and its 
dynamics is usually formulated as an inhomogeneous partial intcgro-differcntial equation. In earlier work [5|-|7[ we 
investigated the relation between this equation of motion and classical phase space trajectories. However, inspired by 
Feynman's Lagrangian description of quantum mechanics Q , in this paper we have explicitly constructed a propagator 
(i.e., a Green's function) for the phase space distribution. 

Consider first a single-particle system, for simplicity in the notations in one dimension, whose coordinates are denoted 
by x. The quantum-mechanical amplitude for the system to go from position x a at time t = t a to position Xb at time 



* Corresponding author: dries.sels@ua.ac.be 



t = t b is given by the Feynman path integral 



x(t b )=x b 



K (x b ,t b \x a ,t a ) = I Vx (t) exp [-S [x (t)} 



(III) 



where S [x (t)] is the action of the system for a trajectory x (t). Let {\n)} denote a complete orthonormal set of states. 
The amplitude A (rib, tb\n a , t a ) for being in state \rib) at time tb, if initially in state \n a ) at time t a , is then given by 



A(n b ,t b \n a ,t a ) = J j (n b \x b ) K (xb,t b \x a ,t a ) (x a \n a ) dx a dx b . (II. 2) 

The corresponding transition probability P (rib, tb\n a , t a ) — \A(rib,tb\n a ,t a )\ 2 can thus be written as 

P (n b ,t b \n a ,t a ) = (x' b \n b ) (n b \x b ) K* (x' b ,t b \x' a ,t a ) K (x b ,t b \x a ,t a ) (x a \n a ) (n a \x' a ) dx a dx b dx' a dx' b . (II. 3) 



Therefore the total probability to be in state \n b ) at t = t b is the sum over all possible transitions from \n a ) to \n b ) , 
weighted by the initial probability to be in a given state \n a ): 



P(n b , t b ) = ^2 p ( n a, ta)P (n b , t b \n a ,t a ) , 



(II.4) 



such that 

P(n b ,t b ) = j j J J (x' b \n b ) (n b \x b ) K* (x' b ,t b \x' a ,t a ) K (x b ,t b \x a ,t a ) |^ P(n a ,t a ) (x a \n a ) (n a |x' }| Ax a Ax b dx' a Ax' b , 

^ (IL5) 
where the term between braces is identified as the initial density matrix p (x a , x' a ,t a ) of the system, since 

p (x, x', t)=^2 (x\n) P(n, t) (n\x') . (II.6) 



Because of the orthonormality of the states, the probability to be in a state \n b ) at time tb can also be extracted from 
the density matrix by 

P(nb, h) = J J (x' b \n b ) (n b \x b ) p (x b , x' b , t b ) dx b dx' b . (II. 7) 

It follows, by comparison with Eq. (|II.5l) . that the density matrix at time t b is given by 

p(x b ,x' b ,t b ) = J J K (x b ,t b \x a ,t a ) K* (x' b ,t b \x' a ,t a ) p(x a ,x' a ,t a )dx a dx' a . (II. 8) 

The density matrices can be expressed in terms of the Wigner distribution function by an inverse Weyl transform 

f(x, P ,t)= j P L+i, x -i, t y-rt/K^L o P (x,x',t)= i f(^±^,p it y+^-z')^dp. (n.9) 

Defining the propagator K w of the Wigner function as 

f(x b ,Pb,tb)= J J K w (xb,Pb,tb\x a ,Pa,t a ) f (x a , Pa, t a ) dx a dp a , (11.10) 

one may extract it from Eq. (|II.8[) : 



K w (xb,p b ,t b \x a ,p a ,t a ) 



/// 



x(t b )=x b ri(t b )=£ b 



(t a )=x a J$(t a )=t 



VxVte*v[- [ -Kit* +8 



x4 

X+ 2 



-S 



X -~2 



(11.11) 

2-Kh 

(11.12) 
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Note that the explicit time dependence in the path variables x (t) and £ (t) is omitted in order not to overload the 
notations. 

For example, consider x to be the coordinate of a (nonrelativistic) particle with mass m in an external potential 
V(x), described by the action 



S[x]= I ( — x 2 - V{x) ) dt. (11.13) 



One then finds 



mx^-V [ a;+| ) + V (x- f I J dt. 



2/ V 2 , 

After an integration by parts for the kinetic term one obtains the following expression for the propagator K u 

rx(t b )=x b r£(tb)=£b 

>x{t a )=x a h(t*)=a* 



rr rx{t b )=x b r£(tb)=£ b 
K w (x b ,p b ,t b \x a ,pa,t a ) = // / / DxT)^ exp I — (mx — p) 

J J Jx(t a )=x a Je<t a )=f a V" 



X {~h L \r* +V T + 2 )- V \ X -V) At J 2*h 

The path integral over £ sums over all possible paths between two fixed points (£ aj £b) , but additionally one has to 
integrate over all possible initial and final points (£ a , . Performing this last integral results in an unconstrained path 
integral over all possible paths. Moreover, the integral over the boundary term fixes the initial and final momentum 
p to mi, such that one is left with 

x(t b )=x b . 

K w (x b , Pb ,t b \xa,Pa,t a ) = ^ J** J™** Vx J V£cxp (~ J ' (mx(;+V (x+?\ -vU-^j dt) . (11.14) 



P(ta) 



-mx,, 



Note that in previous work we obtained the same expession (Eq. (7) in Q) for the propagator of the Wigner function 
by an infinitesimal time lapse expansion of the Wigner-Liouville equation. With this alternative derivation we have 
shown that the Wigner function can also be constructed directly from the Lagrangian form of quantum mechanics 
without the need of the usual canonical quantization procedure. This is important because some problems are much 
more complicated in Hamiltonian form, or a Hamiltonian description might not even exist. Furthermore, we have 
shown that the propagator (|II.14|) for the Wigner function is completely general, and not restricted to the phase space 
evolution of pure states. 

It is clear that the double path integral 14|) can only be solved analytically for a limited set of problems. As an 
example of such a problem consider a particle in a time dependent quadratic potential V(x, t) = a(t) + b(t)x + c(t)x 2 . 
In that case the propagator can be explicitly evaluated: 



K w (x bi p bl t b \x a ,Pa,ta) = I ,. w T>X \ £>£ 
quad J X{ta)—Xa J 



p(£& f cx ? (-1 It! (™+ b (*) + 2c (*) x) fdt) 



p{t a )—mx a 



x(t b )=x b 

oc f p[ * h)=mXi Vx S (mx+b (t) + 2c (t) x) 

J q(t a )—x a 
p(t a )—mx a 

= S(x b - X cl (t b \x a ,p a ,t a )) S (p b -p c l (t b \x a ,p a ,t a ))), (11.15) 

where x c \ (t b \x a ,p a , t a ) and p c \ {t b \x ai p a ,t a ) are the position and momentum at time t b along the classical trajectory 
with initial position x a and momentum p a at time t a . This confirms the well known result 0, 0) that the Wigner- 
Liouville equation can exactly be solved by the method of characteristics for all harmonic problems. It is the basic 
formula for most of the subsequent calculations. In the appendix a detailed derivation is given, based on the Feynman 
propagator for this type of potentials. 



III. INFLUENCE FUNCTIONALS IN PHASE SPACE 



In the previous section we derived the phase space propagator for a single-particle system. But many interesting 
systems consist of interacting subsystems. One of those subsystems is usually of particular interest. In this section 
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we generalize the previous phase space description with the help of influence functionals Q , allowing to describe the 
behavior of the subsystem of interest, coupled with an external (quantum) system, solely in terms of its own variables. 

Consider a more complicated system that consists of two subsystems with coordinates x and u. The subsystems are 
coupled by a potential Vi(x,u), which is incorporated in the action as Si [x,u] . 

The extension of the path integral thus takes the form 

i>x(t b )=x b pu(t b )=U b / ^ \ 

K (x b ,u b ,t b \x a ,u a ,t a ) = / / VxVuexp I - (S x [x] + S u [u] + Si [x, u]) ) (HIT) 

Jx{t a )=x a Ju(t a )=u a \ n / 

where S x [x] and S u [u] are the non-interacting contributions of the subsystems x and u to the action. Again, for 
brevity, the path variables x and u are implicitly assumed to be time dependent. If {\n)} denotes a complete 
orthonormal set of states for the x subsystem, and {\j)} similarly for the u subsystem, the amplitude (|II.2j) can be 
generalized to 

A(n b ,j b ,t b \n a ,j a ,t a ) = ffff( nb \ Xb } Ub\ u b) K (x b ,u bl t b \x a ,u a ,t a ) (u a \j a ) (x a \n a ) dx a dx b du a du b , (HI. 2) 

The conditional transition probability P (n b , j b ,t b \n a , j a ,t a ) = \A (n b , j b , t b \n a , j a , t a )\ 2 for subsystem x to go from 
state \n a ) at t a to state \n b ) at t b , while subsystem u goes from \j a ) to \j b ) thus becomes 



P{n b ,jb,tb\n a ,ja,t a ) = JJJJ JJJJ (x' b \n b ) (u' b \j b ) K* (x' b ,u' b ,t b \x' a ,u' a ,t a ) (j a \u' a ) (n a \x' a ) x 

x {n b \x b ) (j b \u b ) K (x b , u b , tb\x a ,u a ,t a ) (u a \j a ) (x a \n a ) dx a dx b du a du b dx' a dx' b du' a du b . (III. 3) 

From hereof we assume that only the description of subsytem x is physically relevant. In other words, one is interested 
in the probability P (n b , t b \n ai t a ) of the transition of subsystem x from state \n a ) to \n b ) . This can be found from the 
conditional probability (|III.3[) by summing over all final states \j b ) and initial states \j a ) , weighted by the probability 
that subsystem u was initially in state \j a ) ■ 

P (n b ,t b \n a ,t a ) = }~] P(j a ,t a )P(n b , j b , tb \n a ,j a ,t a ). (HI. 4) 

ja ,jb 

Substitution of the conditional probability (|III.3[) into Eq. (|III.4|) and regrouping terms leads to 

P (n b ,t b \n a ,t a ) = \\\\ JJJJ ( x 'b\ n b) K* (x b ,u' b ,t b \x' a ,u' a ,t a ) (n b \x b ) (n a \x' a ) K (x b ,u b ,t b \x a ,u a ,t a ) (x a \n a ) 

^2 (u a \j a ) P(j a ,t a ) (ja\u' a ) > < ^ (u' b \j b ) (j b \u b ) > dx a dx b du a du b dx' a dx' b du' a du b . (III. 5) 



3a I I Jb 



The sum over the initial states \ j a ) is clearly identified as the initial density matrix of subsystem u, whereas the closure 
relation ensures that the sum over all final states \j b ) reduces to 8(u b — u' b ). Then, substitution of expression (|III.1[) for 
the propagators K (• • • ) and rearranging terms one ends up with the following expression for the required transition 
probability: 

P (n b , t b \n a ,t a ) = JJJJ ( x 'b\ n b) (n b \x b ) (n a \x' a ) (x a \n a ) x 

rx(t b )=x b rx'{t b )=x' b / ■ \ 

x / / ■DxVx'exp(-(S x [x}-S x [x'}))J r [x,x'}dx a dx b dx' a dx' b , (111.6) 

Jx{t a )=x a Jx'(t a )=x' a \fi / 

with 

r u{t b )=u b 



F[x,x'}= JJJJ p(u a ,u' a ,t a )S(u b ~ u' b ) \J Vuexp ^ (S u [it] + S r [x,u]) 

C ru' (t h )=u' b / ■ \ 1 

x I J Vu' exp y— (S u [u'\ + Si [x', u'])J \ du a du b du' a du' b . (111.7) 
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J- [x, x'\ contains a double path integral as indicated by the braces, and it can be regarded as an influence functional Q 
since it describes the full influence of subsystem u on subsystem x. 

The relation (|II.9|) readily allows to write the influence functional in terms of the initial Wigner distribution function 
of the u system. Expressed in the center-of-mass and relative coordinate system it becomes: 



T [x, x'] = 



u(t b ) = U b /•r]( t b)=Vb 



T>uDrj exp 



i / S u [u + %] + Sj 



S u [ 



X, M+f ] 



iSV \x',u— 51 



2ttH 



du a dr/ a du b dr] b , (Uli 



where a factor <5(?7b) was replaced by its plane wave representation. 

Proceeding along the lines (|II.3I III. 8|) as in section [TT] one obtains the time evolution of the reduced density matrix 
of subsystem x: 



p(x b ,x' b ,t b ) = 



x(t b )=x b fX '(t b )=x[ 



x(t a )—x a Jx'{t a )—x' a 



VxVx' exp ( - (S x [x] - S x [x']) ) T [x, x'] } p (x a , x' a , t a ) dx a dx' al (III.9) 



which differs from (|II.8|) merely by the occurrence of the influence functional in the path integral. 

Note that we have assumed the two subsystems to be initially independent so that the probability P{j a ,t a ) of 
finding u in state \j a ) is independent of the state of x. This means that the initial total density matrix was supposed 
to be separable. 

Some useful properties of the influence functional arc listed below. The identity 



F* \x,x'} = T[x',: 



(111.10) 



follows directly from Eq. (IIII.7|) by interchanging x and x' . Note that this property will later ensure that the propagator 
for the Wigner function is always a real quantity. Furthermore, if there are a number of statistically and dynamically 
independent subsystems uj acting on x and if ^(x^x') is the influence functional of the jth subsystem on x, then 
the total influence function is the product of all the individual functionals : 



N 



J r [x,x']=Y[j 7j [x,x'\ 



(III.ll) 



This property is a direct consequence of the total initial density emerging as a simple product when all uj are 
statistically independent. If they are also dynamically independent, then each density matrix can propagate separately. 
Finally, it is often convenient to write the influence functional in the form 



T [x, x'] = exp — $ [x, x'] 
n 



(111.12) 



where $ \x,x'] is called the influence phase. For independent subsystems, as in 1 1|) . the corresponding influence 
phases add. According to Eq. (|III.10|) . the influence phase turns out to be antisymmetric under the exchange of x 
and x' if the phase is real, and symmetric if the phase is imaginary. More properties on influence functionals can be 
found in @. 

It follows from Eq. (|III.9|) . and in analogy with the analysis presented in section [TT1 until Eq. (|II.12[) . that the 
propagator for the reduced Wigner distribution function becomes 



K w (x b ,p b ,t b \x a ,Pa,t a ) 



X(t a )=x a J£(t a )=e 



VxV^ exp 



/ 


( 






\ 










i 




+s x 


x+l 


T 


x- 






h 




















V 


V 


-s x 


x — | 


) 


J 











If the action of the system x is of the form 



d&djq 

(111.13) 



(y i 2 ~ V(xj) dt, 



(III. 14) 
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(i.e., without a magnetic field) the further analysis in section [TTI until Eq. (III.14[) simplifies this propagator: 



K w (Xb,Pb,tb\x a ,Pa,t a ) = 7TT 



x{t b )=x b 
p(t b )=mx b , 



2ttH J x(t a )=x a 
p(t a )—mx a 



Vx / ViT 



x+ r x ~2 



xexp,-- 



X ~l 



at 



(111.15) 



The crucial ingredient for further development is the influence functional J- [x, x'\ . Clearly representing the propagator 
of subsystem u under the influence of subsystem x, the path integrals in (|III.7|) can be calculated analytically if 
S u M + Si [x,u] is quadratic in u. Below the results for a bare harmonic action S u with a linear coupling Si will be 
discussed in some detail. 



IV. HARMONIC SUBSYSTEMS WITH LINEAR COUPLING 
A. A single oscillator 

In this section we have reduced the subsystem it to a single harmonic oscillator with mass M and with potential 
energy V (u) = Mu> 2 u 2 /2, interacting with the subsystem x of interest. Taking the interaction energy to be u ■ j(x), 
with an arbitrary function j(x ) and a linear dependence on u, we are left with 5/ [x,u] = — f* b uj t {x)dt, where we 
added a subscript t to 7 to remember the time at which its path variable x should be evaluated. The total action of 
subsystem u, to be used in the influence functional (|III.8[) . is thus 

tb / M .0 M 



S u [u] + Si[x,u] = J \—u-—u z u 2 -u lt {x)'j&t. (IV.l) 
The argument in the exponent of Eq. (|III.8|) thus becomes linear in 77 and in u: 

T[x,x']= J J jjjjf^P^) ei{Pa ^ PbVb)/h j^ b _ ^ Vue^^- l -j\{ lt {x)- lt {x'))&t^ x 

r M=Vb „ ( * f tb f M .. M 2 it(x) + j t (x') \^\ d Pa d Pb j 

x J Vrjcxp y-- J \Murj - Moj z ur)-r) Jdtj ^ du a dr/ a du b drj b . 

After an integration by parts of the kinetic term J* b iifjdt = Ub"i]b ~ u a rj a — f* b ilr]dt, and imposing Mu a fi = p a .b, the 
path integral over r\ becomes unconstrained: 

u(t b )=u b / _ 



T\x,x'\ = JJJJ f(u a ,Pa,t a ) J^J^Vuexp (^\j t {lt{x)- lt {x'))udt*j X 



p(t a )=Mu a 

x / P^exp (i jf* (-Mil M^u- lt{x ' )+ 2 lt{x) ^ rjd?j d P« d P^ du ». (iv.2) 

The path integral over all r\ in the last line is of the form of the quadratic path integral pi.!5[) . and restricts the phase 
space trajectories of u 

M(u + u;\)+ lt(x)+ 2 lt(x,) =0, (IV.3) 

with the formal solution 

\ u a . \ /"' ls(x) + l s {x') sinw (t — s) , 
u{t)=u a cosLu{t-t a ) + —smLo(t-t a )~ / '" K ' n ' sy ' i J -ds. 

The remaining path integral on the first line of PV.2[) imposes that the initial velocity is u a = p a /M. The trajectories 
in u are thus reduced to a single path, with the conditions p a b = Mii a & at the end points, which also eliminate the 
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integrations over Ub and pb- One thus readily arrives at 



t h rt 



F[x,x'\ = cxp y 2hMuj J J (ls(x) +j s (x'))(j t (x) -j t (x'))smoj(t- s) dsdtj x 

f (u a ,p a ,t a )exp (uacosuj (t - t a ) + -j?— sin u> (t - t a )j (7t(x) -7t(x'))dtJ dr a dp a7 (IV.4) 

One might be concerned about the normalizing factors accompanying the several delta functions in the derivation, 
but this possible problem is resolved by considering the uncoupled limit 7 = 0. An alternative derivation, using the 
explicit solution of the Feynman path integral 0] for the action (S u [u] + Si [x, u]) , confirms this result. 

Clearly the first line in expression PV.4I) is independent of the initial state of the harmonic oscillator. This term 
represents an effective interaction of system x with itself. The second line is an expectation value which transfers all 
necessary information about the initial state of the harmonic u system into x. 



1. Example: initial wave packet 

Despite the classical trajectories which govern its dynamics, the influence functional PV.4P is intrinsically of quan- 
tum mechanical nature, because the initial Wigner distribution function / (u a ,p a ,t a ) of the oscillator is bound to 
satisfy the uncertainty principle. A sharply defined initial distribution like S(u a — UQ)5(p a — po) can not be of the 
form pi. 91) . However, a valid initial wave function could be a Gaussian wave packet 

G( , a) v / a7^ V 4A2 / ' 

From the corresponding density matrix pa {u,v! ,t) = (u', t) ^>g (u, t) one easily finds the Wigner distribution 
function pi.9[) of this wave packet 

fc(u,p,t a ) = -e W ^ — 2A j. (IV.5) 

Then, by performing the integrations in Eq. pV.4p and rearranging some terms one finds the following influence 
phase piI.12j) : 

$G [x, x'} = — j^rj j j ( ls (x)+ ls {x , )){ lt {x)- lt (x'))sinu{t- s) dsdt 

~ u oJ (lt{x) - j t (x)) cos uj(t-t a )dt- -j^j- I (j t (x) - j t (x )) sinu (t - t a ) dt 



(A 2 + . a 'ifrl i I COSCJ (s — t) 

( 7 ,(x) - %(x')) ht(x) - lt {x')) A 4A h2 M \J dsdt (IV.6) 

^Ju Ju \ + (A 2 - 4A 4 V cos a; (s + 1 - 2t a ) 



The real part of the influence phase is the same as one would obtain from a (forbidden) initial Wigner distribution 
function S(u a — uo)S(p a — Po)- Besides the effective interacting of system x with itself, the real part now contains an 
external driving potential which oscillates in time with a frequency ui and its spatial dependence is given by j(x). The 
magnitude of this driving potential depends on the initial average (vacuum) displacement of the oscillator (uo,Po) ■ 
The imaginary part of the influence phase results from the uncertainty on the initial position and momentum of u; 
its physical significance becomes more apparent in the next example. Note that some terms vanish under specific 
conditions, e.g., if the average initial position or momentum are zero. The last term becomes zero whenever A 2 = jfitui- 
This condition is satisfied if the oscillator was initially in an unsqueezed coherent state. In that case one finds the 
ground state or vacuum influence phase: 

$ va c [x, x] = J J (%(x) + 7 S (x')) (7t(x) - 7t(x'))sinw (t - s) dsdt 

t b 

( 7s (x) - j s (x')) ( 7f (x) - lt {x')) cosw (s - t) dsdt. (IV. .7) 



AMuj 
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2. Example: thermal equilibrium 

The case of thermal equilibrium at the start of course deserves some additional attention. The initial equilibrium 
Wigner function Q of the harmonic oscillator u is then given by 

tanh hBtuj ( tanh \Bfojj ( „ „ n 2 \ \ 

f eq (Ua,Pa,t a ) = exp f j^— ( Mw 2 u\ + J^)l (IV.8) 

It follows from Eq. pV.4|) that the equilibrium influence phase $ e<? [x, x'] is 



2Mu 



? e 9 [x,A = — i— / / (j s (x)+ ls (x')){j t (x)-j t {x'))smoj(t-s)dsdt+ 



icoth-Bhcu f tb f tb 

(lt(x) - lt {x')) (-f s (x) - 7s (x'))cosw(i - s )dsdt (IV.9) 



AMuj 



ta J t 



2 ) ~~ 



While the real parts of $ e? and $„ ac are the same, the imaginary part of $ e? has increased by a factor coth 

as compared to <fr va c, &s a consequence of thermal broadening of the distribution. Note now that coth 

[1 + 7is(w)] + nsipS), with rig (a;) the Bosc-Einstcin distribution with zero chemical potential. From a physics point 
of view one would therefore associate the imaginary part of $> eq with the emission and absorption of the quanta 
represented by the oscillator. When the temperature is zero system x can only interact with the zero point fluctuations 
of u, allowing only losses trough spontaneous emission. A more detailed discussion on this can be found in 

B. Many independent oscillators 

The generalization of the result from one to many oscillators is trivial if the oscillators are independent. According 
to property (jIII.ll|) the total influence functional then becomes the product of all the individual influence functionals. 
Consider every oscillator Uj to have a mass Mj, a frequency u)j and consider the interaction energy with x to be 
Uj ■ 7j(x). Then we immediately arrive at the following expression for the total influence phase of a collection of N 
oscillators 

F N [x,x']=exp^^2 2M uj, / / ^ j <yXs ^ + 13 ^ ^ <Cl3 ^* ^ ~ lj ^ ^ ^ sin ^ ^ ~ ^ dsd ^j X 

N r r /„•/■* 

X 



\\ jjfi (Ua,Pa,t a )cxp jf ^ a COSWj (t - t a ) + SmUJ 3 ( f _ *»)^ (lj( X t) ~ Jj( X 't)) d ^ dr a dp Q , 

(rv.io) 



3=0 



where fj (u a ,p a , t a ) represents the initial Wigner function of the jth oscillator. Let us consider a simple example now 
that can easily be generalized to more complicated situations. 

1. Example: Equilibrium oscillators with bilinear coupling 

If we assume the interaction energy to be bilinear in {v,j, x} , such that u ■ Jj(x) = ^jU ■ x, and if additionally all 
oscillators are initially in thermal equilibrium pV.8[) . then the influence phase associated with influence functional 
(IIV.IOP is given by 

N „,2 ,t„ ft 



—j — / / (x s + x' s ) (x t - x' t ) sinwj (t ~ s) dsdt 
j=0 AMjUj J ta J tn 



4M,-w 

3=0 J 



3 — cotYi-fiftkjj I I (x t — x' t )) (x s — x' s )) cos ujj (t — s) dsdt. 



ta J t 
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In the continuum limit, when TV — > oo while uij+i — tOj —> 0, we can assume there is a distribution of oscillators, such 
that the relevant weight Y(lj)&uj of the oscillators between to and to+dto is 

it 



3=0 



In this case the influence phase becomes 

rib rt / poo 



$eq,many[x,x'] = -J^ (J T(u) SmU (t - s) dijj (x s + x' s ) (x t ~ x' t ) dsdt 

tb r t b r poo 



fhT(to) , Bhw 
— — coth - — 
o V 4 2 



cos to (t — s) dto 



(x t - x' t )) (x s - x' s )) dsdt. 



If x is the coordinate of an otherwise free particle with mass m wc find the following propagator for its reduced Wigner 
function 



K w {Xb,Pb,tb\x a ,Pa,t a ) = 7T-r 



x(t b )=x b 
1 f p(tb)=mx b , 



2ttH J x(t a )=x a 
p(t a )—mx a 



R(t - s)£ t ■ Zsdsdt 



x cxp 



A(t - s)x s ds ■ £ t di 



with A(t) and R(t) defined as: 



A(t) = / T(io) sin(w£)dw, 



R(t) 



^Mcothf^ 
4 V. 2 



cos(tot)duj . 



2. Example: Caldeira-Legett model and thermalization 

In the specific case that T(to) = 2rjio, as considered in detail by Caldeira and Leggett in we obtain A(t) = 
—2r]8'(t). If additionally the temperature Tb is high enough such that we can approximate \^y~ coth f^rp^ ~ 
4 + 0(/3), wc find i?(t) = r)kTbS(t — s). Therefore we arrivcQ at the following expression for the propagator in the high 
temperature limit 

x(t b )=x b t \ / ■ pt 

K w (xb,Pb,tb\xa,,Pa,t a ) — — |^ Jl^ =x ' J exp (--^- J £ t 2 dt J exp ( - ^ / (mi + 77a?) • &di J . 

p{t a )—mx a * a * a 

It can easily be shown, for example by timcslicing the £ path integral, that this propagator can be rewritten as a 
single path integral: 

x{t b )=x h . 

K w (xb,Pb,tb\x a ,p a ,t a ) = Jx(t]=x b ® xex P f - ^ \.rp b J (mx + r]x) 2 dt 

p(t a )—rnx a ^ a 

If we are not interested in the real space motion of x, for example because the initial distribution of the particle is 
homogeneous in space, but only in the marginal propagator K w (pb,tb\p ai t a ) to go from p a to pb, then we get 

rp(t b )=Pb ( 1 f tb ( (p) 2 m\ 2 P 2 V \ \ 



1 under some debatebly assumptions [3, [Toll about the limits in the integral. 
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An integration by parts shows that the term in pp only contributes at the boundaries: 

The remaining path integral is a simple Gaussian Feynman path integral. This path integral can be solved with 
standard techniques which yield the following expression for the reduced momentum space Wigner propagator: 



K w (p b , tb\p a , to) = , exp 



1 {Pb-Pa CXp (tb ~ tg))) 2 



^TrmJfeT^l-expC-^^-O)) V 2mfcT& 1 - exp (- ^ (t 6 - <„)) ^ 

The maximal transition probability is, as expected, attained along the solution of the classical equation of motion 
p + —p = 0, i.e. when p b = p a exp (— ^ (t b — t a )j . The variance of the propagator is given by 

cr 2 = mkT b ( 1 - cxp ( — — (tb - t a ) 



The time evolution of the effective temperature T e of the particle x is therefore given by 

T e (t b - t a ) = T b (l- exp ( -— (t b - t a 



such that the system thermalizes at a characteristic time r = . Finally consider the time tb — t a 3> t, such that the 
system is thermalizcd and T e ~ T b , then the propagator becomes 

1 Zip 2 
lim K w (p b ,t b \p a ,t a ) = -===exp (--r^- — 

This means that the reduced Wigner function of the particle thermalizes into a Maxwcll-Boltzmann distribution, 
regardless its initial Wigner function. 



V. CONCLUSION 



In conclusion we have used Feynman's Lagrangian description of quantum mechanics to express the propagator of the 
Wigner function as a path integral. Propagating Wigner functions rather than wave functions has a double advantage. 
First of all one can directly propagate uncertain initial configurations in time instead of pure states only. Secondly the 
attained Wigner function propagator becomes a delta function in the classical trajectory for all harmonic problems 
and subproblems.In contrast to our previous derivation in [?J based on the Wigner-Liouville equation, the present 
treatment does not rely on a canonical quantization procedure of the Hamiltonian. This Lagrangian formulation, with 
the help of influence functionals, allows for a transparent description of two interacting subsystems, and to find the 
reduced Wigner function propagator of one of the two subsystems. In the last section we have generalized this result 
to a system interacting with many other, mutually independent, subsystems. As an example we considered in more 
detail the model by Caldeira and Legett, and showed how to calculate the reduced Wigner function propagator for 
it, using the techniques explained in this manuscript. 
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Appendix A: Wigner propagator for quadratic potentials 

In this appendix the propagator 15|) of the Wigner function for quadratic potentials of the form 

V quad (x, *) = o(t) + b(t)x + c(t)x 2 (A.l) 
is derived from the well known [ll| Feynman propagator -fTquad {x b , tb\x a , t a ) for this type of potentials: 



/ 77Z I % \ 

-ft^quad (xb : tb\x a ,t a ) — \ ~ 7Z „ / 7T 6Xp — Squad (x b , tb \x a , to) I , 

y Zmnj [tb, t a ) \n J 



(A.2) 



S qua d (a;&,ii,|a; ,i ) = i 2 - a (t) - 6 (i) .t - c (t) a; 2 ) dt, (A. 3) 

where Squad (a^6, tb\x a , t a ) is the action of the system along a classical trajectory from (x a ,t a ) to (x b ,t b ) , to be 
determined from the classical equation of motion 

mx + 2c (t) a; = -b (t) with * = Xa ' (A.4) 
w w x {t b ) = x b . v 

The function / (tb,t a ) only depends on the initial and final time, and is independent of the positions and momenta. 
It is the solution of the differential equation 

(m^ + 2c (i)) / (t, t a ) = with = i. (A.5) 

For general time dependence of c (£) , these differential equations rarely have a solution in closed form, but it is sure 
that two linearly independent solutions, say x\ it) and X2 (t) , of the homogeneous equations exist: 

mil, 2 + 2c (t) n i2 = 0. (A. 6) 

Their Wronskiaan X2 (t) — dz j t %i (t) is independent of t, because (|A.6|) reveals that its time derivative is zero 

%^x 2 (t) - ^^X! (i) - VF independent of t. (A.7) 

Civ Cvl/ 

Since / (t,t a ) is also a solution of the homogeneous differential equation (|A.6[) . it is a linear combination of X\ (t) and 
x 2 (t) . Taking the boundary conditions into account it becomes 

/ (t, t a ) = k{t ^ a) with h (a, t) = x x (s) x 2 (t) - x 2 (s) x x it) . (A.8) 

If one imposes that the solution of the homogeneous differential equation (|A.6[) exhausts the boundary conditions 
at t a and tb, the trajectory x (t) is of the form 

x = ±Mi + ^2±x b + x p (t) , (A.9) 
h(t a ,tb) h(t a ,t b ) 

where the particular solution x p (t) has to satisfy the boundary conditions x p (t a ) = = x p (tb) . It is easily found by 
the variation of parameters method, with the result: 

PW h(t a ,t b )J ta m W h(t a ,t b )J t m W K ' 

It is fairly easy to calculate the initial and final velocities i a ,b, which are of particular relevance below: 

1 ( dh(t a ,t b ) [*» b(s) \ 

x a ^- x b W + / h(s,t b )ds), (A. 11) 



h(t a ,t b ) \ dt a 

*>> = 7777~TT ( X " W + x b dk{ t tb) - C — h ( 
h(t a ,t b ) V ot b J ta m 
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Applying an integration by parts j x 2 dt — xx — J xxdt in the kinetic contribution to the classical action, it can be 
rewritten as 



Squad (Xb,h\x a ,t a ) = ^ ([xx} t t J i b - [ ^-x(t)dt 



2 



2 dh(t a ,t b ) _ 2 dh{t a ,t b ) 
-b dt b X a dt a 



tb b{s) f tb b(s 



+ Wx a x b -x b h (t a , s) ds — x a / h(s,tb) ds 



n{t a ,tb) \ 2 J ta m J ta m 

b(s)x p (s)ds, (A.13) 



~ 2 

where the boundary velocities and the homogenous contribution x (t) have been filled out. 

In the propagator (|II.12[) for harmonic interactions of the form (|A.1[) . only terms linear in £ Q ^ survive in the 
exponent, and one is left with 

mW 

K w {Xb,Pb,ti)\x a ,p a ,t a ) = - "2 - - X 

quad (27m) h(tb,t a ) 

* h (s* (-+ wht (— t - jC ds ))) * 

x /cxp(i t ( P „-^(,„^|^ -^ + /''^(,Wd»))) dfc. (A.14) 

where the remaining integrals are S functions. Taking the results for the boundary velocities into account, one thus 
readily finds 

K w (x b ,p b ,t b \x a ,Pa,ta) = 5 (x b - X c \ (t b \x a ,p a ,t a )) S (p b - Pd {t b \x a ,p a ,t a ))) , (A. 15) 

quad 

which is the desired result 15|) . 
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